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Abstract In Nachtergaele obtained explicit lower bounds for the spectral gap 

f""^ ■ above many frustration free quantum spin chains by using the "martingale method" . 

We present simple improvements to his main bounds which allow one to obtain a 
sharp lower bound for the spectral gap above the spin-1/2 ferromagnetic XXZ chain. 
As an illustration of the method, we also calculate a lower bound for the spectral gap 
of the AKLT model, which is about 1/3 the size of the expected gap. 
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^ ' 1. Introduction 



In this paper, we obtain lower bounds on the spectral gap above frustration free 
ground states, which are special ground states which exist for certain special quantum 
i -^h | spin systems. In the simplest "frustration free" quantum spin chain can be 

defined as follows: For each site, let Tt x = C d be a finite dimensional Hilbert space, 
such that all 7i x are equivalent to one another. For any finite subset A C Z, one 
defines Tt\ = ® x( z\7~(- x . Consider an interaction $, defined on all finite subsets of 
Z, such that ^(Ao) £ B(Ha ) f° r au finite Ao. Then for any finite subset A C Z, we 
^ ■ define 

H(A) = $ ( A o) > 

AoCA 

where, as usual, ^(Ao) £ B(H\ a ) is extended to Ha by tensoring with the identity 
on Ha\a . This is the set-up for a typical quantum spin system. What distinguishes 
a frustration free quantum spin system are two assumptions 

(1) For each Ao, ^(Ao) is a positive operator. 

(2) For each nonempty A, the subspace £(A) = Ker(iJ(A)) satisfies dim((?(A)) > 
0. 

In this case, we call Q (A) the ground state space, and declare the ground states to be 
frustration free. 

We point out an obvious implication. Since all ^(Ao) are positive, 
Ker(JT(A)) = n AoC A Kcr($(A„)) . 
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This means that any ground state ip £ Q(A) not only minimizes the total energy, i.e. 
is a minimum energy eigenvector of -ff(A), it also minimizes all the local interaction 
energies. This is to be contrasted with a frustrated spin system. 

To mention just a few interesting examples of frustration free systems, there are 
ferromagnets: the isotropic Heisenberg ferromagnet has a frustration free ground 
state, as do the anisotropic "kink" ferromagnets known as XXZ models. A frustration 
free antiferromagnet is the AKLT model, invented by Affleck, Kennedy, Lieb and 
Tasaki [jlj , which is an important model for understanding the role of spin-dimension 
in antiferromagnetic chains. (Haldane ]l2| conjectured that the antiferromagnetic, 
isotropic Heisenberg chain has gapless excitations above the ground state if the spin 
is half- integer but is gapfull if the spin takes an integer value.) There are numerous 
stochastic models, such that the Markov generators can be viewed as frustration free 
models. Of these, we just mention one, the Kac model for which the gap has 
recently been computed exactly by Carlen, Carvalho and Loss B using a technique, 
similar in spirit, but different than the martingale method, which is specific to mean 
field models. We point this out to indicate that this is still a very active field of 
research. 

Henceforth we will assume that $ is only a nearest neighbor interaction, and that 
it is translation invariant. Neither of these requirements is essential.^] This allows us 
to write Ha = J2 X x+ieA h( x > x + 1); where h(x, x + 1) is a fixed positive operator of 
just two spins. 

2. Main Theorem 

The main result of this paper is a method of calculating lower bounds for the 
spectral gap, by diagonalizing a family of finite dimensional matrices, with dimensions 
much smaller than the starting Hamiltonians H\, A C Z. There are two situations 
where this is particularly useful: First, if there is a symmetry which allows explicit 
diagonalization of these matrices, which is the case for the spin-1/2 XXZ model. 
Second, for the finitely correlated states (see Section |^) the matrices involved are all 
of a fixed finite dimension. If the dimension is small enough, the matrices can be 
diagonalized by hand, which is the case for the AKLT model. 

Suppose A = [a, 6] and Ao = [c, d] C A. Then we define H{c, d) — Ha , interpreted 
as an operator on Ti.\. We suppress the dependence on A from the notation. Similarly 
we define G(c, d) = Ker H(c, d). We define G(c, d) — Proj(C/(c, d)) to be the orthogonal 
projection onto Q(c, d). The spectral gap, 7(0, d), is defined to be the largest number 
such that 

H(c,d) > 7(c,d)(l - G(c,d)) . 
This number is independent of A, and depends only on c — d + 1 = #[c, d\. We 
define 7/v = min2<„<jv 7(1, n), which has the useful property of being a nonincreasing 
function of N. A sensible convention is that H(a,a) = 0, Q(a,a) = H a = C d , and 
G(a,a) = 1. However, there is no way to define 7(0,0), since both H(a,a) and 
1 — G(a, a) vanish. 

The auxiliary Hamiltonian, whose analysis is easier than that of Ha and which gives 
information on 7(0,6), is just the ground state projector restricted to an appropriate 

2 An example of a non-translation invariant frustration free model was treated in [^j. There, the 
XXZ model is perturbed by a defect, which nevertheless is still frustration free. In that paper, the 
"martingale method" is applied in a slightly different way because there the spin chain is grown on 
both edges, while we will grow the spin chain along one edge only. 
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subspace. Suppose Ao = [— to, ti] C A= [a, b]. Define 

(2.1) e(m,n) = sup{ (ip\G(Q,n)ip) : ip € Q(-m,0) , ip S G^m.n) 1 - , ||^|| = 1}. 

The value of e(m, n) is independent of a and b as long as a < —to and b > n, and 
by translation symmetry the answer would be the same if the interval [— to, n] were 
replaced by [x — m, x + n] for any x. For m,n > 1, we define 

£m,n = sup e(m', n) , 

m'~>m 

which is monotone decreasing in to, but not necessarily in n. 
Our main result is the following: 

Theorem 2.1. Suppose that s m ^ n < 1/2, for some n > to > 1. For any N >2 and 

ip G ^(l,^)- 1 , de/me 

ip! = (1 - G(1,N - n))1> , ^ 2 = G(1,N - n)ip , 
whence it follows ip = ipi + ip2- (If N < n + 1 t/ien V^i — 0, t/>2 = V'J TTien 

(2.2) 0>\H(1, N)i>) > 7 m+ n(a(£ m: „)||^i|| 2 + P{e m , n )U 2 \\ 2 ) , 



where a(e) = (yl — e — and /3(e) = yl — i(yl — £ — y^)- 

Our improved bound is due to the following simple estimate. 

Lemma 2.2. (a) Suppose H is a positive matrix, G = Proj(Ker(iJ)), and 7 is the 
gap; so that H > 7(1 — G). If (j) is a vector which is not in the kernel of H - not 
necessarily orthogonal to the kernel either - then for any ip 

(2.3) (VW) > - G)4>)\ 2 / {<P\{1 - G)4>) . 

(b) If G is an orthogonal projection, <p is a nonzero vector, and ip is any vector 
orthogonal to (p, then 

(2-4) |(^|G0)| 2 <|H| 2 ||0|| 2 (G) (1-(G),), 

where (G) = (0|G0)/||0|| 2 . 

Proof: Both are by Cauchy-Schwarz. First, 

|(^|(1 - G)0)| 2 = |((1 - G)m - G)0)| 2 < ||(1 - G)^|| 2 ||(l - G)^|| 2 

<7- 1 (^|^)ll(l-G)0|| 2 . 

For the second, choose an orthonormal basis <pi, 4>2, ■ ■ ■ with <p\ — <p/\\<f>\\. Then 

1/2 / \ 1/2 



MGct>)\ < ]T M<M(<i>n\G<t>)\ < (ElW^I 2 ) (^r(G</#„)(^|G0> 

n>2 n>2 n>2 

= ||^|(^<G<^ n )W>n|G0) Xl 2 



Moreover, 

^(G^|0„)(0„|G0) = (G<f>\G<i>) - {GflfaXfalGt) = ||G0|| 2 (1 - ||G0|| 2 /||0|| 2 ) . 



n>2 



Proof: (Of Theorem 2.1) The proof is by induction. The first step is to check that 
the proposition is trivially satisfied for N < m + n. Indeed, since a < [3 < 1, we have 

> max(a7 m+ „, f3j m +n) , 



4 



by the definition of "f m +n- The next step is the induction step. 

For the induction hypothesis we suppose the proposition is true whenever TV < No, 
where No > m + n, and then prove the proposition is also true for N = Nq + 1. This 
means we must show 

(2.5) (^\H{l,N)^) > lm+n [a\\M 2 + PU2\\ 2 ] 

where G(l, N)tp = 0, Vi = [1 - G(l, N - n)]ip and ip 2 = [G(l,N - n) - G(l,N)]ij>. It 
is useful to make a further orthogonal decomposition 

ibn = [1 - G(l,N - 2n)]Vi , V12 = [G(1,N - 2n) - G(1,N - n)]Vi . 

(We have assumed that N > 2n, otherwise we set ipu = 0.) By the induction 
hypothesis, we conclude 

(^|JT(l,JV-n)Vi) > 7m+»H| Villi 2 +^ll^i2|| 2 ]- 

Since tp 2 G Q{1,N - n), it is clear that (ip\H(l, N — n)ip) = (<0i|if(l,iV- n)0i). So 
to prove equation (|2.5|), we just need to prove 

(2.6) (^ji?(7V - n, N)1>) > 7m+n(/3||0 2 || 2 - (fi - a)\\i> 12 \\ 2 ) . 
We assume 

(2-7) ||<A 2 || 2 > (l-^)||0 12 || 2 = -^L||<M 2 



because otherwise equation (2J3) is automatically satisfied. (Henceforth, we will write 
e in place e m . n .) 

To estimate the left hand side of equation (2.6), we use Lemma 2.4. Specifically, 
by part (a), substituting ip 2 for <j>, we have 



(iP\H(N -n,N)^} > 7m+n J 



[1 - G(N — n, N)]ip; 



2 



2 



(i} 2 \[l-G(N-n,N)]ii> 2 ) ' 
where we have used "f m +n < 7n+i- We define the number 

r, = (G(iV - n, N))^ = (^ 2 |G(7V - n, N)^ 2 )/U 2 \\ 2 
which allows us to rewrite 

(2.8) (iP\H(N - n,N)i>) > 7m+»l(#l - G(iV - n, N)}^ 2 )\ 2 / (1 - r))\\^ 2 



The denominator is not zero because, by (2/7), ipi ^ 0, and by the definition of £ m ,n, 
V < £m,n < 1/2. 

Writing G in place of G(N -n,N), 

(2.9) (0|[1 - G]0 2 ) = (1 - 77) || V^2|| 2 - (<M|GV> 2 ) - (0 12 |G^ 2 ) . 
It is easy to see that 

(2.10) (ipn\G^ 2 ) = (ij n \G(N-n,N)ip 2 ) = 

because ipn and G(N — n, N)ip 2 are both eigenvectors of the Hermitian operator 
G(1,N — 2n), with different eigenvalue. (Note that G(N — n, N) commutes with 
G(l, N ~ 2n) because they are localized on disjoint intervals.) By Lemma 2.4 (b), 



(2.11) |(0i 2 |GV> 2 }| < HV12II • \\HWv(l-ri)- 
Combining equations ( |2.9| ), ( 2.1 0| ) and ( 2.11 ), gives 

(2.12) (0|[i - G]^ 2 ) > V^v\\M(V^v\\M - VvWM 



Using and observing that < J^_ e , the right hand side of ( 2.12 ) is positive. 
Hence 



(2.13) |<#1 - G}^)\ 2 > (1 - V^fiV^vWHl - 0?ll^i2|l) 2 ■ 

Therefore, by (|2.8|), equation (pl3) would follow from 



(2.14) (Vl^ll^ll - vWdl) 2 > PU2W 2 -(J3- a)\\i> 12 \\ 2 . 
The proof is reduced to checking that the bilinear form 

B(x, y) = {fl-a + v)* 2 ~ 2^(1 - rj)xy + (1 - 77 - (3)y 2 
is positive. The conditions for positivity are 

(2.15) /3-a + ?7>0, 1 — 77 — /3 > , 

(2.16) 77(1 — 77) < (1-77-/3)03- a + r/). 

They must be satisfied for all 77 in the range [0, s] (where we write e for e m , ra ). Con- 
dition ( 2.16|) is equivalent to (1 — 0)(f3 — a) > r)(2f3 — a). By condition ( fc.lEp , this is 



the same as 

*?< (1 -/?)(/? -a)/(2/3- a) 
for all r\ < e. However the inequality with 77 = e subsumes all those with r\ < e. Again 



using (2. IE), this is equivalent to 

a </3(l- 2e- f3)/{l-e- P). 
The maximum value of a satisfying these constraints occurs when 



a= (VI - e - Vef , P = 1 - e - y/e(l - e) , 
which is just what we set out to prove. | 

The previous theorem can be used to obtain uniform bounds for the spectral gap 
above all finite volume ground states. The situation for infinite volume ground states 
is usually more complicated than just taking the limit of the finite volume gaps. An 
infinite volume ground state is a state lo on the quasilocal observable algebra Ao 
which is locally stable in the sense that for any finite volume A, and any X G Aa, 
u>(X*S(X)) > 0, where S is the derivation S(X) — \im\> [Ha> , X] . For such a ground 
state, one can define the gap 7(0;) to be the greatest number 7 satisfying 

uj(6(X)*S 2 {X)) > yu(5{X)*5(X)) , 

for any local observable X £ A\, with A finite. The purpose of taking three powers of 
5 on the left hand side is to restrict to the positive energy states, since many ground 
states may be quasilocal perturbations of one another. With this definition, we can 
state the main corollary. 

Corollary 2.3. Suppose that e m ,n < 1/2, for some n > m > 1. Then there is a 
uniform lower bound on the finite volume spectral gaps 



inf 7(JV) > 7 m +n(l ~ 2y / e m ,„(l - e m . n )) 



Moreover, if u) is an infinite volume ground state, obtained as a weak-* limit of finite 
volume ground states, then j(u>) obeys the same bound. 
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Proof: The proof of the first part follows trivially from Theorem 2.1. One simply 



observes that (3(e) > a(e) and a(e) = 1 — 2ye(l — s). 

The proof of the second part follows an argument of Koma and Nachtergaele [[L7| , 
which was stated for the XXZ model. Suppose lu is a weak-* limit of finite volume 
ground states lon(' ■ ■ ) = (V'ivl" • • V'Jv)) where the unit vectors ipN are each in Ti\ N 
for An / Z. Then for any local X £ ,4a, there exists an No with X £ A\ N for 
all N > iVo, and wpf) = limjv^oo cjat(A), the limit starting at iV = No- Now note 
that S(X)*S 2 (X) is also local, localized in A^ 2 ) consisting of all points no more than 
2 units from A. There is also an No such that Ajv D A^ 2 ) for all N > No, and then 
we evaluate 

u N (S(X)*6 2 (X)) = (8(X)iI> n \H(A n )5(X)iI>n) . 
Let ijj' N = S(X)ipN. Then if (j> is any ground state of H = H(A^), 

(4>\il/ N ) = (4>\(HX — XH)i/jn) = 
because Hipff = Hcf> = 0. So ip' N is orthogonal to the ground state space, hence 
uj N (5(xyS 2 (X)) > j(N) W N \i// N ) = j(N) lu n (S(X)*S(X)) . 

Using the uniform lower bound for the finite volume gap, and taking liminf's proves 
the claim. ■ 



Theorem 2.1 apparently is a stronger result than the last corollary. For exam- 
ple, one can use it to prove that, whenever the lower bound on the infinite volume 
gap is achieved, then the first excitation is not an "edge state". The only example 
we know where the lower bound is attained is the XXZ model, and there the first 
excitation is actually not an eigenvalue but the bottom of a continuous band of spec- 
trum with quasi-momentum 0. This is an important difference between the XXZ 
model and the AKLT model, where the lowest variational excitations (see [[[o)) have 
quasi-momentum equal to n. 



3. Example 1 : The XXZ model 

The XXZ ferromagnet, with fcmfc-inducing boundary fields is a translation invari- 
ant, nearest-neighbor model with anisotropic exchange 
h(n, n + 1) = - sech(£)S„ ■ S n+1 - (1 - sech(£))S^ +1 + ■ tanh ( e )( 5 3 _ ^ 

where £ > determines the amount of anisotropy. The isotropic Heisenberg model is 
£ = 0, while £ = oo is the Ising model with kink boundary conditions. For j = 1/2, the 
Hamilton is invariant under the action of the quantum group SU 9 (2), where q = e~^. 
In this section we will prove the following proposition. 

Proposition 3.1. Let £ > 0. Then the spectral gap above any infinite volume, weak-* 
limit of finite volume kink ground states is 7(£) = 1 — sech£. 

The representation of SU g (2) is generated by three spin operators 

N N N 

S 3 (l,tf) = £ 5 n' 5+(l,iV) = £e-(»- 1 )^+, S-(l,N) = J2e- iN - n)( S-. 

n— 1 n— 1 n— 1 

The first is the usual generator of U(l); the other two are deformations of the spin 
raising and lowering operators for SU(2). Notably, they are not permutation invariant. 
For us, the fundamental feature is a decomposition property. Define S 3:+, ~(a,b) = 
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r a 1 (5 3,+ ' (1,6 + 1 — a)), where r is the translation automorphism. Then for any 
1 < n < N 

S + (l,N) = S+(l,n) + e- n ^S + (n + l,N), 

S~{1, N) = e- {N - n ^S-(l, n) + S~(n + 1, N) . 

The representation theory 8115(2) for q > is identical to that of SU(2) in many 
respects. In particular, the Clebsch-Gordon series is the same. This means the irre- 
ducible representations of each dimension, or multiplets, occur with the usual multi- 
plicities in H(1,N). Particularly, there is a unique maximal dimensional multiplet. 
As in the isotropic case, it is known [^2[ |[ [llj that the ground state subspace of the 
XXZ fcrromagnet is precisely this maximal dimensional multiplet. Moreover, it was 
shown that all the infinite volume ground states are obtained as weak-* limits 



of these states, which means Proposition 3.1 is particularly relevant. We refer to a 



generic n-dimensional SU g (2) multiplet as V(n). 

To prove the proposition, we use Corollary [^3[ Nachtergaele already calculated 
every e(m, 1) in p0| . We will repeat his argument, for completeness. This also 
serves as a slightly simpler preparation for the calculation in the next section. A 
major motivation for this letter is the fact that Nachtergaele's calculation, with our 
improved bound, is sufficient to calculate the exact gap. 

We now state a lemma which is generally helpful for calculating e(m,n). 

Lemma 3.2. The Hermitian operator K(—m,n) = G{— m, 0)G(0, n)G{~ m, 0) has 
range contained in the Hilbert space Q(—m, 0) <8>£?(1, n). The eigenspace for eigenvalue 
1 is precisely £7(— m,n). The number e(rn,n) is the largest eigenvalue less than 1. 

Proof: The first statement is true because G(0, n)G(—m, 0) has range contained 
in £7(1, ro), and [G(— m, 0), G(l, n)] = implies K(—m,n) also has range contained in 
£7(1, n). Clearly £7(— m, n) is an eigenvalue-1 eigenspace of K(—m, n). If K{—m,n)i^ — 
i/j, then G(— to, 0)ip — ip and (V>|G(0, n)ip) = 1, which implies G(0,n)ip — ip. Since 
if is a nearest- neighbor Hamiltonian, this implies ip € Q(—m,n). The last statement 
follows from the minimax principle. | 



Proof: (Of Proposition 3T) In the case of the XXZ model, the space Q{— to, 0) C 
m, 0) is the unique (m + 2)-dimensional multiplet V(m + 2). The space £7(1, 1) 
is just V(2). By Clebsch-Gordon, V(m + 2) <g) V{2) = V(m + 3) V(m + 1). 
Since G(— to, 0) and G(0, 1) are SU 9 (2) invariant, so is K{— to, 1). This implies that 
K(—m, 1) = Ai Proj(V r (m + 3)) + A 2 Proj(y(m + 1)). Of course, since V(m + 3) is 
the maximal dimensional irreducible multiplet in H(—m, 1), we know Ai = 1, by the 
lemma. We can restrict to the two-dimensional (invariant) subspace of Q{— m,0) ® 
£7(1,1), which is the intersection with the eigenspace {ip : S 3 (—ni,l)t[j — m/2-0}. 
Then it is easy to determine A2 from the trace. In the following we will use a 
nonorthogonal basis, which simplifies the algebra a bit. 
Every vector can be written as 

= (aS-(-m,0) + /3e- (m+1)€ S-(l,l)|fi) 

for some a, f3, where |f2) is the all up-spins state. Note 5 ,_ (1, 1) is the spin lowering 
operator for the interval [1, 1] = {1}; i.e., it is the spin lowering operator at the site 
1. We observe that 

rj> = (aS-(-m, -1) + a e - m< ^-(0, 0) + /3e~ (m+1) « S~ (1, 1)|0) . 
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Since S~(0,1) = S~(0,0) + e-^S~(l, 1), it is clear that 



G(0, 1)V = (aS-(-m, -1) + ^^""^-(0,0) + e -^-(l, 1))) |fi) 
Similarly, since S~(-m, 0) = S"(-m, -1) + e _mi S r_ (0, 0), we have 

G(-m, 0)G(0, 1)^ = f (Z + e- 2 "*)- 1 ^ + e ~ 2 ^ ^± e 1^ )g-(- m , 0) 
V 1 + e ^ 



e ^. e -(m+^s-{l,l))\n) 



1 + e" 2 « 

where Z = \\S-{-m, -1)|^)|| 2 = (1 - e - 2m «)/(l - e - 2 «). (We are exploiting the 
fact that G(0, 1) and G(—m, 0) are rank-one projections restricted to the appropriate 
subspace.) Hence the eigenvalues are equal to the eigenvalues of the matrix 



1 



Z(l + e-^) + e- 2r < e -2(m+i)€ 
(Z + e- 2m «) e- 2 «(Z + e- 2m «) 



(Z + e- 2m «)(l + e- 2 «) 

The trace of this matrix is 1 + e~ 2 «Z/((Z + e~ 2m «)(l + e~ 2 «)). Since one of the 
eigenvalues is 1, that means the other eigenvalue is 

e~ 2 «Z e- 2 « 1 - e- 2m « 



e(m, 1) = 



(Z + e- 2 '< ) (1 + e- 2 « ) 1 + e- 2 « 1 - e -2(m+i)€ ' 

This is increasing with m. So, £1,1 = e~ 2 ^/(l + e~ 2 ^) < 1/2 for £ > 0. Also, 72 = 1. 
So 72(1 - 2 v /e M (l-E M )) = 1 - sech£. 



By Corollary 2.3 this gives the lower bound for the gap of Proposition 3.1. But it 
is also easy to obtain this as an upper bound for the gap above the all-up-spins state, 
using a sequence of variational wave function approximating a long spin- wave state of 
quasimomentum zero. Since the all-up-spin state is a weak-* limit of the kink ground 
states upon translation to the left, this upper bound also works for the gap above the 
kink states. ■ 

Note: The exact gap of the spin-1/2 XXZ chain was previously determined rigor- 
ously, by different means, by Koma and Nachtergaele . 

4. Example 2 : The AKLT model 

In this section we will show how to calculate e(m,n) for the AKLT model. The 
AKLT model is a spin-1, translation invariant, nearest neighbor model with interac- 
tion 

h(n, n + 1) = P 2 (S n + S n+ i) 
where Pi is projection onto the spin-2 states in the tensor product of the two spins-1. 
This model was introduced by Affleck, Kennedy, Lieb and Tasaki, in p|] . They calcu- 
lated the ground states, proved that the ground states are exponentially clustering, 
and proved the existence of a positive spectral gap [Q. However, their existence proof 
gave no explicit lower bound on the gap, which is is a question of physical interest 
(see e.g. In |pl]| , Nachtergaele developed a version of the martingale method for 
quantum spin systems to obtain explicit lower bounds on the gap of many finitely 
correlated spin chains. Finitely correlated spin chains are a class of frustration free 
quantum spin chains, including the AKLT model [^J. The martingale method was 
originally developed by Lu and Yau to obtain the exact scaling of the spectral 
gap and logarithmic Sobolev inequality constant for certain stochastic spin systems. 
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The key point of Lu and Yau's proof for the spectral gap was the introduction of a 
filtration, such that the Gibbs measure is a martingale with respect to that filtration. 
This is the same idea for the case of quantum spin systems, except in place of the 
Gibbs measure one uses the ground state. The martingale property guarantees that 
the decomposition one obtains for the Hilbert space is actually a direct sum decompo- 
sition. In retrospect, this same decompositon is also the key element of the existence 
proof in Q, although they did not call it a martingale decomposition. 

There are several versions of the frustration free ground states of the AKLT model. 
Affleck, Kennedy, Lieb and Tasaki first wrote them down in the Ising basis, as a 
weighted sum over certain configurations which are "alternating sign rows" . These 
are configurations of +'s, — 's and O's such that between any two +'s there is a — 
and between any two — 's there is a +. Aside from this rule, O's may occur anywhere. 
Physicists refer to this as diluted Neel order. On the basis of this, it was conjectured 
(7) that the AKLT state has a long range order, for an order parameter which itself 
is a product of arbitrarily many terms. Partly in order to explain this, Kennedy and 
Tasaki [|l4| came up with a nonlocal unitary transformation. One of the nice features 
of their transformation is the fact that it makes the ground states into product states. 
This simplifies calculations, so we shall use their representation. 

There is another representation of the AKLT ground state, which we shall not use, 
but we would like to mention, which is the quantum Markov chain (known as matrix 
product representation by physicists). Fannes, Nachtergaele and Werner || used this 
to prove many interesting features of the valence bond solid states, in particular that 
all the correlations could be modeled on a finite-dimensional space. They called states 
with this property, which include many generalizations of the original VBS states, 
finitely correlated states. More recently, physicists Kolezhuk and Mikeska Jl6j have 
used this approach to study spin ladders, and derived new examples of frustration 
free models beyond those which were known for spin chains. 

The Kennedy- Tasaki unitary transformation is U = Yij<k exp (z7r S'JS'^). It makes 
the four ground states of the AKLT model, which are otherwise somewhat complicated 
states, into product states <E"fc(l, N) = where <j>\p = [|0) ± \/2\+)\/yJZ 

and 563,4 = [|0) ± V2|-)]/V3. We denote $i(a,b) = ® b n=a (</>,:)„. These four states 
are not orthogonal, but we can calculate the Gram matrix 



M(a,b) := [($,(«, 6)|$>,6))]t, =1 = K^) 6 "^ 1 ]^! 



Specifically, 

(4.17) M(l, N) = h + 3 _JV 



(-1)* 1 1 
{-l) N 1 1 

1 1 {-l) N 
1 1 {-1) N 



We denote the inverse as W(a, b) := M(a, b) 1 , and it equals 



W(1,N) = (l + 2(-l) iV 3 



No-N 



(4.18) 



L-3 



-N 



3-3 

-(-1) N 2 
(-1)" 

1 

1 



-2jV\-l 



-(-1)*2 
1 
1 



1 
1 

-{-l) N 2 
(-1)" 



1 

-(-1) N 2 
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We can write the ground state projector as 

4 

G(a,b) = Ma,b)W(a,b)^(a,b)^ 

where $>j(a,by is the covector dual to $j(a, b). Indeed this operator is clearly Her- 
mitian, and it has the correct range. It is an easy exercise for the reader to check that 
it behaves correctly on the four ground states. 

We will use Lemma ^2 to calculate s(m, n). An important point is that we do not 
need to use an orthogonal basis to obtain a matrix representation. So, instead, we 
use the basis 

0y (— m, n) = m, 0) <x> %(l,n) . 

We will calculate the matrix for G(—m, 0)G(0, n)G(—m, 0) in this basis. Of course 
G{— to,0) fixes @ij(~m 7 n). Before calculating G(0, n)6y(— to, n), we note that for 
any a < c < b, $i(a, b) = $i(a, c) ® $ 4 (c + 1, b). Hence, 

G(0, n)Qij (—to, n) 

= J2 n ) ® (^C . »)«*i(0. ^) t $ l (-TO, -1) <g> $,(0) <g> ^(1, n)) 
= ]T P7(0, n) w M(0, 0) K M (1, n) w #i(-m, -1) <8> $ fc (0, n) . 

kl 

Finally, 

G(-m, 0)G(0, 7i)e y (-to, n) 

= E E **-(-m, 0) <g> (w(-m, 0) rs 17(0, n) w M(0, 0) K M (1, n)y 

rs fcZ 

$ s (-to, 0)t$ i (-m, -1) ® $ fc (0, 0) <8> $ fc (l, n)" 
t7(0, n) w M (0, 0) M M(1, n)y W(-m, 0) rs M(0, 0). fc M(-m, -l) si e rfc . 



E 

klrs 



Let ^4 be the 16 x 16 matrix with 

A(ij, rk) = £ W(0, n) kl M (0, 0) H M (1, n) y W(-m, 0) rs Af (0, 0) sk M (-m, . 

is 

The top four eigenvalues are 1, corresponding to the four- fold degenerate ground state 
space Q(—m,n). Therefore, e{m,n) is the fifth largest eigenvalue. 

To diagonalize A, we observe that while U has obscured the SU(2) symmetry of 
Hn-i it is still there. It's just no longer a "local" symmetry. Each of Q(— m, 0) and 
G(l,n) is equal to 7(3)© 7(1). Thus, the space on which A acts is (7(3)©7(1))® 2 = 
7(5)©37(3)©27(1), using Clebsch-Gordon. Abusing notation, we refer to 37(3) and 
2F(1) as the spin-1 and spin-0 subspace. A is SU(2) symmetric since each G(a,b) is. 
Using this and Schur's lemma, we know that 7(5) is an eigenspace of A with eigenvalue 
A5 G Q. We also know that the eigenspace of A for eigenvalue 1 is Q(—m,n) = 
7(3) © 7(1). Let 7(1)' be the orthogonal complement of G(-m, n) n 27(1) in 27(1). 
Since this is one dimensional, it is an eigenspace for A with eigenvalue Ai € Q. Let 
W be the orthogonal complement of G{—m, n) n 37(3) in 37(3). Then A restricted 
to this block, and restricted to an eigenspace of S 3 , is a 2 x 2 matrix. There are two 
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eigenvalues A3,±, which are irrational, but solvable by radicals. So A is solvable by 
hand, and the eigenvalues are 

1 (l-(-l/3) m )(l-(-l/3)") 



A 5 (m,n) 
Ai(m, n) 
A 3± (m,n) 



9 (1 - (-l/3) m + 1 )(l - ' 

1 (l-(-l/3) m - 1 )(l-(-l/3)"- 1 ) 

9 ' (1 - (-l/3) m + 1 )(l - (-1/3)"+!) ' 

1 1 + (-1/3)" 1 + (-1/3)" - 19(-l/3) m +"+ 2 ± 2y/R(m, n) 



9 (l-(-l/3)™+ 1 )(l-(-l/3)"+ 1 ) 
R(m,n) = [(-l/3) m - (-1/3)™] 2 - 128(-l/3) 2m+2n+4 

+ 4(-l/3) m+n+2 (l + (-l/3) m )(l + (-1/3)™) . 
To calculate £i,i, we set n — 1 and observe that the maximum eigenvalue occurs for 



m 



= 1, A 5 (l, 1) = 1/4. Since 7(2) = 1, Corollary 2.3 gives 



inf j(N) > 1 - « 0.133975 . 

N>2 ' y ' ~ 2 

If we instead use £3,3, we must numerically diagonalize the six-site chain to calculate 
76- We have done this using Lanczos, and the answer is 76 = 0.3985(±0. 00005). 
Taking n = 3 we find e(3, 3) = A 5 (3, 3) = (7/20) 2 . So we obtain 

inf 7(JV) > 0.137194. 

AT>2 

Several authors have calculated variational upper bounds for the spectral gap above 
the AKLT ground states. One of these authors was Knabe Jl5| who not only calculated 
variational upper bounds, but obtained rigorous lower bounds, too. Fath and Solyom 
|l0[ showed Knabe's trial states to be equivalent to the single mode approximation 
states of Arovas, Auerbach and Haldane and themselves gave a third interpretation 
as hidden domain walls. These all give an upper bound of j(N) < 10/27 = 0.370 .... 
Our lower bound is not spectacular, but at least has the correct order. By yet different 
methods, Fannes, Nachtergaele and Werner obtained the excellent lower bound 
y(N) > 3/10. 
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